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1. Introduction 


The 2001/77/CE European Commission Directive sets the target 
of 22% of gross electricity generation from renewables for the 
Europe, by 2010. In the case of Portugal the target is more ambitious 
at 39%. In order to satisfy this target, it is planned the installation of 
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additional generation capacity from renewable sources, in which the 
largest part of the increase it is in wind power, whose target is to 
reach 5700 MW by 2012, and 8000 MW by 2020. 

In a scenario of large scale penetration of renewable production 
from wind and other intermittent resources, it is fundamental that 
the electric system have appropriate means to compensate the 
effects of the variability and randomness of the wind power 
availability. Portugal has a relatively large hydro capacity of 
4600 MW and only about 50% of the potential is used. Large solar 
photovoltaic plants have been installed and are planned. The 
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foreseen sharply decreasing costs will allow a large increase of the 
solar capacity in the next 10-15 years. To ensure an optimal mix in 
the medium-term, three main renewable sources (hydro, wind and 
solar) deserve to be considered. Also, through the proper 
application of demand-side management (DSM) and demand 
response (DR) technologies it is possible to influence the load 
shape to mitigate the unbalance between intermittent power and 
the demand. 


2. Wind power intermittence 


Wind energy has characteristics that differ from conventional 
energy sources. If the contribution of this production vector in 
energy terms is not a cause of concern, the contribution to the 
power balance, and therefore the impact in the supply security 
needs attention due to the intermittent and random character of 
this production option. 

The output of wind power is driven by environmental 
conditions outside the control of the generators or the system 
operators. Since the wind is determined by random meteorological 
processes it is inherently variable. Supply of power from wind 
turbines is stochastic in nature and the actual power is more or less 
proportional to the third power of the wind velocity (from cut-in 
wind speed to the rated wind speed). The wind output varies 
seasonally between summer and winter [1] and the variations are 
also present on shorter time scales, namely on hourly basis [2]. 
Unlike conventional capacity, wind-generated electricity cannot 
be reliably dispatched or perfectly forecasted, and exhibits 
significant temporal variability. 

Several extreme ramp rates were recorded during storms [3]: 


e Denmark—2000 MW (83% of capacity) decrease in 6 h or 12 MW 
(0.5% of capacity) in a minute on 8th January 2005. 

e North Germany—over 4000 MW (58% of capacity) decrease 
within 10 h, extreme negative ramp rate of 16 MW/min (0.2% of 
capacity) on 24th December 2004. 

e Jreland—63 MW in 15 min (approx 12% of capacity at the time), 
144 MW in 1h (approx 29% of capacity) and 338 MW in 12h 
(approx 68% of capacity). 

e Portugal—700 MW (60% of capacity) decrease in 8 h on 1st June 
2006. 

e Spain—800 MW (7%) increase in 45 min (ramp rate of 1067 MW/ 
h, 9% of capacity), and 1000 MW (9%) decrease in 1 hand 45 min 
(ramp rate —570 MW/h, 5% of capacity). Generated wind power 
between 25 and 8375 MW have occurred (0.2% and 72% of 
capacity) in a single year. 

e Texas, USA—loss of 1550 MW of wind capacity at the rate of 
approximately 600 MW/h over a 2.5 h period on 24th February 
2007. 


In addition to being variable, wind power production is also a 
challenge to accurately predict on the time scales of interest to: day 
ahead and for long-term planning of system adequacy. It is possible 
to forecast in a concrete zone, the average wind power density for 
the whole year; however, it is impossible to precisely forecast the 
days or the hours with wind [4]. 

Integrating large amounts of wind energy into the electric 
generation mix requires some special considerations. Beyond the 
variability, a lot of wind generation occurs in hours when energy 
use is low. The uncontrollable nature of wind makes it less valuable 
to system operators than dispatchable power. The variability and 
uncertainty of wind energy production require that power system 
operators take measures to manage its delivery. These measures 
may increase the cost incurred to balance the system and maintain 
reliability. Substantial amounts of wind generation in a utility 
system can increase the demand for the various ancillary services. 


Previous results also reveal a diminishing benefit as wind power 
penetration increases [5]. 

In 2005, the impact of meeting 100% of Western Denmark’s 
annual electrical energy requirement from wind energy was 
analyzed [6]. The results of the study demonstrated that the 
system could absorb about 30% energy from wind without any 
wasted production. However when wind share reaches 50% the 
excess wind energy starts to grow considerably. With the 
hypothesis of total energy demand of 26 terawatt-hours (TWh) 
generated by wind power, 8 TWh of the wind generation would be 
surplus because it would be generated during periods of low 
consumption. In that case, the electricity costs would double, 
because of the need to use other costlier back-up systems. 


3. Options for managing intermittency 


The connection of wind turbines to the electricity grid can 
potentially affect supply reliability and power quality, due to the 
unpredictable fluctuations in wind power output [7]. In a scenario 
of large scale penetration of renewable production from wind and 
other intermittent resources, it is essential that the electric system 
have appropriate means to compensate the effects of the 
variability and randomness of the renewable power availability. 
This concern was traditionally addressed by the promotion of the 
wind resource studies and by the identification of solutions based 
on reversible hydropower dams [8]. However, in the electric 
system planning, other options deserve to be evaluated. 

The intermittency of wind energy can be reduced by some 
techniques: 


e Grid integration. 

e Technical distribution of the generators. 

e Geographic distribution of the generators. 
e Improved forecasting techniques. 


These techniques have as aim the increasing of the predictabil- 
ity of the production and the substantial reduction of the global 
variations. However, although those improvements bring benefits, 
several periods of low wind production and substantial variations 
will remain. Thus tools to respond to short- to medium-term and 
long-term variability will be necessary, managing the operational 
and capacity reserve, respectively. For large scale integration of 
wind power the provision of flexible capacity reserve will be of 
crucial importance. To achieve that aim several options are 
possible: 


e Power plants providing operational and capacity reserve. 
e Interconnection with other grid systems. 

e Curtailment of intermittent technology. 

e Dispatchable distributed generation. 

e Energy storage. 

e Use of complementarily between renewable sources. 

e Demand-side management. 

e Demand-side response. 


All these options have as aim continuously balancing demand 
and supply and backing up other capacity shortage. The present 
paper analyses the application of the three last methodologies in 
Portugal. 


4. Complementarity between renewable sources 


To evaluate the intermittence and complementarity from 
renewable energy production, a climate data time series over 50 
years was obtained. The collected information includes the global 
solar radiation (monthly average), the wind velocity (monthly 
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average) and the monthly water inflow in dams. The locations for 
data collection were selected based on the approximation between 
the collected data and the annual variation of the wind power, 
solar photovoltaic and hydropower in all the country. 


4.1. Multi-variable climate model 


With the collected data, a mathematical model was developed, 
to generate random years, enabling the study of the sources 
complementarity for a large number of years. The variable with 
lesser deviation from a regular pattern (daily and seasonal cycles) 
is the solar radiation and thus was tacked as the base variable for 
the model. As it can be observed in Fig. 1, the collected data has a 
distribution close to a normal distribution, and thus the solar 
radiation was generated randomly with the Box-Miller transfor- 
mation (1) 


Smx = [sy + Osx v —2 ln Uj sin(27U2) (1) 


where: Smx is the solar radiation of the model (month x); jus, the 
average value of the solar radiation (month x); Osx the standard 
deviation of the solar radiation (month x); U; and U2 is the 
independent values randomly generated by a uniform distribution, 
between 0 and 1. 

The wind velocity is more uncertain than the solar radiation, 
but more predictable than the water inflow, and thus was then 
considered as secondary variable in the model. However, wind 
speed cannot be determined purely by a random distribution, 
because the two variables are not independent. But since a perfect 
correlation does not exist, it is incorrect to determine the wind 
velocity only with the correlation with solar radiation. Thus, wind 
velocity was determined with two components: random compo- 
nent and component obtained with the wind-solar correlation. 

The wind speed probability density function follows a Weibull 
distribution, but the monthly values variation (from year to year) 
also follows a pattern similar to a normal distribution (Fig. 2). 
Therefore the first component was randomly generated by the 
same process that was used to the solar radiation (2). 


Wrx = Uwr + Owe V—2 In Uj sin(27U2) (2) 


where: W,x is the wind velocity (month X) randomly generated; 
Hwx the average value of the wind velocity (month x); ows is the 
standard deviation of the wind velocity (month x). 

The second component was obtained with the correlation with 
the randomly generated solar energy. In that case two possibilities 
exist: when the variation is larger than the obtained with a unity 
correlation (3) or when the variation is lower (4). 


(Smx ae Sym) (2 = Pxs. w) (Wx = Wyxm) 


Wesx = (3) 
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Fig. 1. Distribution of the, average daily, solar radiation, in January. 
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Fig. 2. Distribution of the wind velocity, in January. 


(Smx m Sxm) Pxs,w (WxMm = Wxm) 
Seu = Sxm 


Wesx = (4) 
where: Wes, is the wind velocity (month x) generated with the 
correlation with solar radiation; pxs,w the correlation between the 
solar radiation and the wind velocity (month x); Wxm, Wxm is the 
maximum and minimum wind velocity (month x); Sxm, Sym is the 
maximum and minimum solar radiation (month x). 

The final value is the addition of the correlated value weighed 
by the correlation and the random value weighed by one minus the 
correlation value (5). 


Winx = Wesx > Pxsw — Wrx(1 — Pysw) (5) 


where Wmx is the wind velocity of the model (month x). 

The water inflow is the most uncertain variable and has a 
different distribution (Fig. 3). Due to such shape, a normal 
distribution generation cannot be used. For each month the data 
was divided in intervals and for each interval a number of random 
values proportional to the probability distribution were generated. 

Due to the larger uncertainty, both the correlations with the 
solar radiation and wind velocity were used. The final value is the 
addition of the correlated values weighed by the correlations and 
the random value weighed by two minus the sum of the 
correlations (6). 


Hmx = Hesx + Pysy + Hewx > Pyw y + Arx(2 — PxsH — Pxw.n) (6) 


where: Hmx is the hydro inflows of the model (month x); Hcsx the 
hydro inflows (month x) generated with the correlation with solar 
radiation; xs, the correlation between the solar radiation and the 
hydro inflows (month x); Hcwx the hydro inflows (month x) 
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Fig. 3. Distribution of the water inflow, in January. 
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Fig. 4. Monthly variation of the solar capacity factor. 
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Fig. 5. Monthly variation of the wind capacity factor. 


generated with the correlation with wind velocity; oxw,4 the 
correlation between the wind velocity and the hydro inflows 
(month x); H,x is the hydro inflow (month x) randomly generated. 


4.2. Climate data analysis 


In a way similar to wind power, hydropower and solar power 
are also intermittent resources, due to the dependence of 
meteorological conditions. However the variations do not occur 
at the same time for the three renewable sources and with a proper 
mix the total output can partially be compensated [9]. 

Using the climate model, a series of 500 years was generated for 
the three variables. The three collected variables have different 
units of measurement, and thus to enable a comparison between 
them, a conversion to an undimensional unit was made, the 
capacity factor. As it can be observed in Fig. 4, the solar radiation 
has small fluctuations relatively to the average year and 
additionally the yearly variation curve does not change in shape. 
For each month the extreme values (Max and Min) are also shown. 

The wind velocity presents high variations relatively to the 
average year, with impact in the yearly variation curve shape 
(Fig. 5). 

The hydro inflow presents huge variations relatively to the 
average year and unpredictability (Fig. 6) preventing a reliable 
estimation. Large storage reservoirs can dampen the fluctuations of 
the water inflows. 

Fig. 7 shows the average monthly capacity factor for each 
variable (wind velocity, solar radiation and water inflow). As it can 
be observed the solar radiation is higher between May and 
September, occurring the opposite with the wind velocity and 
water inflow. Thus, the solar radiation has the maximum value in 
July and the minimum in December. Both the wind velocity and 
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Fig. 6. Monthly variation of the hydro capacity factor. 
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Fig. 7. Yearly variation of the wind, solar and hydro capacity factor. 


water inflow has the maximum value in February, registering the 
minimum in September and August, respectively. 

The wind velocity and the water inflow have average variations 
along the year with a very similar course, having the two curves a 
high correlation (0.98). The solar radiation varies almost inversely 
relatively to the wind velocity and the water flow (correlation of 
—0.7 and —0.66, respectively). That observation indicates that the 
complementarity between solar energy and the pair wind power/ 
hydropower is quite high. Solar energy can then be used to face the 
seasonal variations of wind power. The hydropower is not 
complementary to wind power but due to its similar variations 
is the ideal means to store the excess wind energy to cope with the 
intermittence, using its storage, dispatchable power and dynamic 
response capabilities. 


5. Demand-side management 


Rather than attempting to match power generation to 
consumer demand, the philosophy of load management takes 
action to vary the load to match the power available. Through the 
proper application of demand-side management (DSM) technolo- 
gies [10] it is possible to reduce the need of new installed 
intermittent power to achieve the renewable penetration targets. 

The most critical situations due to high penetration of wind 
power occur in periods with high energy consumption. Thus, the 
DSM technologies have a major role in avoiding critical situations 
due to intermittent power, mainly the technologies with impact in 
the peak load reduction. One serious problem occurs in the 
summer days with high temperatures (elevated consumption of air 
conditioning) with reduced wind velocities. Therefore the energy- 
efficient DSM technologies with high impact in the peak load hours 
are the most important. 
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Fig. 8. Impact of the DSM in the Portuguese load diagram (January 2020). 


Taking as reference the electric energy consumption in 
Portugal, in the year of 2008, an evaluation of the consumption 
evolution was made, considering a business-as-usual scenario 
(BAU) and a scenario with DSM measures application. In such 
evaluation, a consumption increase of 3%/year (BAU scenario) and 
the application of DSM measures corresponding to 1% of the year 
consumption (DSM Scenario) were considered. 

A variety of demand-side management technologies were 
considered to accomplish such objective, trying to achieve a larger 
impact with a minimal cost. The aggregated impact in the load 
diagrams of the selected technologies in the residential (lighting, 
appliances and space conditioning), services (lighting, office 
equipment and space conditioning) and industrial (lighting, energy 
efficient motor systems and drives) sectors was determined. As 
example, Fig. 8 shows the global load diagram (joining the impact 
of the three sectors), on a work day, in January, to 2020. In such 
image the largest impact of the DSM measures in the peak load 
hours can be observed. In this example, a reduction of 8.3% in the 
peak load was obtained, due to the applied DSM measures, 
minimizing the most dangerous situations originated by the 
intermittence of the renewable resources. It must be noted that 
such impact, caused by the DSM measures, is obtained with an 
average cost of 0.023 €/kWh, that is much less than the production 
cost of any renewable energy. 

Other technology that can perform a major role in the integration 
of renewable intermittent power is demand response (DR) [11]. 
With these technologies it is possible to direct or indirectly force a 
consumption reduction in critical situations, in a short time. The idea 
behind DR is that if the marginal peak load price is higher than the 
value that a consumer gets out of the services derived from the 
electricity, he would be willing to modify the demand, if paid the 
peak price or slightly less instead. Traditionally the DR technologies 
were typically used to attend to economical concerns. However, 
nowadays they can be used to improve the system reliability, 
reducing instantaneously the energy consumption to prevent the 
most unbalanced situations, like the problems that result from the 
large space conditioning consumption on days with reduced wind 
velocity. As more customers practice automated price-responsive 
demand or automatically receive and respond to directions to 
increase or decrease their electricity use, system loads will be able to 
respond to, or manage, variability from wind power production. 

Additionally to the DSM scenario, controlling 5% of the peak 
load (nearly 490 MW) with DR technologies it is possible to smooth 
the loads diagram. In Fig. 9 it is possible to see the impact of the DR 
measures added to the DSM scenario. With the aggregated action 
of the measures it will be possible to achieve a peak load reduction, 
in 2020 of 12.9%. 
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Fig. 9. Impact of the DR in the Portuguese load diagram (January 2020). 


6. Renewable mix optimization 
6.1. Intermittence 


The intermittence and complementarity between renewable 
sources can be divided in two different time scales: between years 
or between months in an average year. For the intermittence 
between years, the percentage of the wind, solar and hydro 
capacity factor, which maximizes the global capacity factor in an 
average year, must be determined. Also the difference between 
each year and the average year must be minimized, being such a 
difference less than 15%, since this is the minimum possible value, 
based on a 500 year simulation with the existent data (7). The sum 
of determined percentages must be 100% (8) and each one must be 
positive (9). 


(Cts.c fsa + Ctw.c fwa at Cty.c fpa) 


>0.85 7 
(Cts.C f savg TT Ctw.c fwa ns Cti-C f avg) z l . 
Cts + Ctw + Cty = 1 (8) 
Cts, Ctw, Cty > 0 (9) 


Using the climate model, the conditions are tested for a series of 
500 generated years to maximize the energy generation objective 
(10). 


max(Cts.C f savg + Ctw.c fwavg + Ctu.C f Havg) (10) 


where: Ctw, Cts and Cty is the contribution of the wind, solar and 
hydro component; Cfwavg, Cfsavg and CfHavg the capacity factor in 
average year of the wind, solar or hydro component; Cfwa, Cfsq and 
Cfua is the capacity factor in the year A (1-500) of the wind, solar or 
hydro component. 

The optimal solution for the above problem has a share of: 44% 
in wind, 39% in solar and 17% in hydro production. 

For the intermittence between months, in an average year, the 
percentage of the wind, solar and hydro capacity factor, which 
maximizes the global capacity factor, was determined. Simulta- 
neously, a minimum difference between the global capacity factor 
and the desired variation to achieve the monthly consumption was 
searched. Thus, a condition that such difference cannot exceed 23% 
(again, this is the minimum possible value, based on a simulation 
with the existent data), in any month of the average year, was set 
(11). For such analysis the relation between the consumption in 
each month and the average monthly consumption was consid- 
ered. Also the determined percentage sum must be 100% (12) and 
each variable must be positive (13). In the hydropower case, due to 
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the storage capacity associated with these plants, it was considered 
that the monthly difference relatively to the annual average can be 
reduced by 1/3. 


(Cts.c fsp + Ctw-C fws + CtH.C fp) 


> 0.77 11 
Ctcy(Cts.C f savg F Ctw-C fwavg T Cti-C f Havg) E ( ) 
Cts + Ctw + Cty = 1 (12) 
Cts, Cty, Cty > 0 (13) 


Using the climate model, the conditions are tested to a series of 
500 generated years to maximize the objective (14). 


max(Cts.c f sag + Ctw.c f wang + Cti-C f Havg) (14) 


where: Cfwo, Cfsp and cfrp is the capacity factor in the month b (1- 
12) of the wind, solar or hydro component; Cfwavg, Cfsavg and CfHavg 
the average capacity factor of the 12 months of the wind, solar or 
hydro component; Cte is the percentage of the consumption in the 
month b (1-12). 

The optimal solution for the above problem has a share of: 10% 
in wind, 65% in solar and 25% in hydro production. 

Thus, taking into account the above results, the following 
conditions must be minimized in the expansion of renewable 
capacity, in order to ensure minimum intermittence either 
annually (15) or monthly (16). In such analysis the solar power 
was divided in photovoltaic and concentrated solar power (CSP) 
with storage [12]. 


min[2.56Ppy -= 2.56Pcsp = 1.14Pw = 2.94P4] (15) 


min[1.54Ppy + 1.54Pcsp _ 5Pw _ 2Pu] (16) 


where Ppy, Pcsp, Pw, and Py are the installed power in photovoltaic, 
concentrated solar power, wind and hydro power. 


6.2. Peak load contribution 


In Portugal, taking as a reference the summer and winter peak 
loads in a BAU scenario, the impact of the DSM and DR measures 
was considered. The obtained power was weighed by the desirable 
renewable power energy contribution share and the forecasted 
contribution of each renewable source was subtracted. In such 
objective function the use of biomass power was also considered 
(Pg). Thus, the contribution to the peak load can be optimized by 
(17), to winter and by (18), to the summer. 


min[Ctr(Pgauw — Posmw — Pprw) — CtpvwPpv — CtcspwPcsp 


— CtwwPw — CtuwPy — CtawPs] (17) 


min[Ctr(Pgaus — Ppsms — Pors) — CtpvsPpy — CtcspsPcsp 


CtwsPw — CtysPy — CtgsPa] (18) 


where: Ctr is the renewable power contribution to the global 
energy production; Pgayw, Pgaus, Posmw, Ppsms, Pprw and Ppsr the 
peak load in to winter (W) and summer (S) to the BAU, DSM 
reduction and DR reduction scenarios; Ctpyw, Ctpys, Ctcspw, Ctcsps, 
Ctww, Ctws, Ctuw, Ctys, Ctgw, Ctgs is the contribution to the winter 
(W) and summer (S) peak load for each renewable technology. 

The contribution of each power generation technology to meet 
the peak load was forecasted by the interception between the daily 
capacity factor diagram, for each technology, with the hour of the 
peak load, for the winter and the summer peaks, taking into 
account the technologies dispatchability. 


6.3. Production costs 


To optimize the renewable mix, the production costs need to be 
minimized (19). 


min[CpyPpyhpy + CcspPcspltcsp + CwPwhw + CuPuhy 
+ CgPghg] (19) 


where: Cpy, Ccsp, Cw, Cu, Cg is the levelised cost of energy for each 
renewable energy; hpv, hcsp, hw, hy, hg is the operation hours of 
each renewable energy at full power. 


6.4. Restrictions 


In addition to the previous objectives for the renewable mix 
optimization, several restrictions need to be considered: 


e Energy consumption share to be ensured by renewable sources 
(20). 

e Maximum potential for each renewable technology ((21), as 
example to the wind power). 

e Actual installed power for each renewable technology ((21), as 
example to the wind power). 

e Yearly growth to the installed power for each renewable 
technology ((22), as example to the PV power). 


Ppvhpy + Pcspltcsp + Pwhw + Puhy + Pghg 


= Ctr(Egau — Eps) (20) 
Pwi < Pw < Pw, (21) 
Pwi + CrrPw, — Pw > 0 (22) 


where: Egay is the annual energy consumption for a BAU scenario; 
Egau the annual energy consumption reduction due to the DSM 
measures; Pw; the installed wind power; Pwp the maximum 
potential in wind power; Crp is the yearly maximum growth to the 
installed power in each renewable technology. 


6.5. Global optimization 


Taking as base the available climate data, a climate model was 
constructed to generate random years (Section 4), used to test the 
conditions to minimize the intermittence between years and 
months. Taking as a reference the BAU scenario for the energy 
consumption growth, the impact of DSM and DR measures was 
added (Section 5), reducing the required energy consumption to be 
supplied by renewable sources, the non-guaranteed peak load 
share, and the productions costs. Such conditions will optimize the 
renewable mix, subject to the renewable potential, the initial 
installed renewable power capacity, the maximum renewable 
capacity yearly growth and the target percentage consumption 
that needs to be ensured by renewable energies (Fig. 10). 

The multi-objective model of the renewable mix optimization 
have as optimization functions the Eqs. (15)-(19), and as 
restrictions the Eqs. (20) and (21) (example for wind power) 
and (22) (example for wind power). Using the described model, a 
simulation of the years between 2010 and 2030 (with a 5 years 
interval) was performed for the Portuguese data, considering an 
increasing renewable share in the total energy production, from 
40% in 2010 to 60% in 2030. For the optimization process, a weight 
of 50% for the costs component and 50% for the reliability 
component was used. The reliability component was equally 
divided between the 4 components (intermittence between years, 
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Fig. 10. Multi-objective model. 


2010 2015 2020 2025 2030 
== Hydro eee Biomass 


Fig. 11. Optimized Portuguese renewable mix between 2010 and 2030. 


intermittence between months, contribution to the winter peak 
load and contribution to the summer peak load). 

Fig. 11 shows the results of the multi-objective model. In the 
above simulation time period, the obtained results for the objective 
functions were in the direction of an increasing reliability, due to 
the reduction of the intermittence (between years and months) 
and of the non-guaranteed peak load share (winter and summer). 
In the case of the peak load, the non-guaranteed part was reduced 
to 34%, for the summer and 32% for the winter. Such evolution was 
obtained in first place due to the use of dispatchable technologies 
(biomass, concentrated solar power with storage and hydropower) 
and after that due to the use of technologies with high 


complementarity with the wind power (photovoltaic and concen- 
trated solar power). 


7. Conclusions 


To face the intermittence of renewable sources several options 
can be considered. The option to mix complementary energy 
sources like wind power, solar power and hydropower will 
mitigate the problems, when comparing with only using one 
source of renewable energy and the use of dispatchable renewable 
sources has a great importance. Using historic data for the wind, 
solar and hydro availability a multi-objective model can be used to 
support such analysis, in order to optimize the complementarity 
between renewable sources. 

DSM and DR can also have a major role, either by reducing the 
needs of new intermittent capacity or by adjusting the consump- 
tion in real time, to face production variations. Thus the demand 
consumption reduction due to the application of those technolo- 
gies should be incorporated in the analysis. 

The described multi-objective model enables the optimization 
of the renewable mix, ensuring a minimum level of intermittence 
(between years and months), a minimum non-guaranteed peak 
load share (both in summer and in winter) and a minimum global 
cost, using the complementarity between renewable sources and 
the cost-effective impact of DSM and DR measures. 
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